DOCOHENT RESaHE 



ED 197 992 



SE 03tt 167 



TI TIE 



INSTXTnTTOR 

S^ONS AGENCY 
POB DATE 

ANT 
NOTE 



Cowan r Christina E- 

Pressure aud Bnoyancy in Aquatic Ecosystems. Physical 
Processes in Terrestrial and Aquatic Ecosystems^ 
Transport Processes- 

Washington Oniv. r Seattle* Center for Quantitative 
Science in ForestrVr Fisheries and Wildlife, 
National Science Fcundat lo-^i , Hashingtonr D- C. 
Jlay 79 

NSF-GZ-29B0: NSF-SE07a-1 76^6 

83p^: For related documents^ see SE 03U 160-166 and 
SE*033 5B1-597. Contains marginal legibility in 
coiuDuter print^outs- 



EDPS PBICE 
DESCSIPTOT?S 



IDENTIFTE'^^S 



n^O^/VCOtx Plus Postage. 

Bio-loay : Biophvsics: ^College Science: ^Computer 
Assisted Instruction: Computer Programs : Higher 
Education: *Tnterdiscipiinary Approach; Learning 
Modules; ^Mathematical Applications : Physics : 
* Pressure fPhysics) : Science Education: * Zoology 
* Aquatic organisias 



ABSTRACT 

This module is part of a series designed to be used 
bv life science students for instruction in the application of 
physical theory to ecosvstem operation. Host modules contain computer 
programs which are built around a particular application of a 
physical process* This module explores some of the characteristics of 
aquatic organisms which can be explained in tt-rms of pressure and 
buoyancy. Emphasis is placed on physical l^ws which illustrate how 
organisms interact with their environment. The two physical 
properties of primary concern are density and viscosity. The effects 
or pressure caused by the density o' water and the probieias of 
flotation and buoyancy in static water are discussed. A problem set 
illustrates and extends ideas presented in the text. An accompanying 
computer procrram, BOOY, permits the student to investigate the 
sinkina of plankters. A knowledge of calculur and differential 
equations is presumed. (Author/CS) 
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PREFACE 

This module explores some of the characteristics of aquatic 
organisms which can be explained in terms of pressure and buoyancy. 
Emphasis is placed on physical laws which illustrate how organisms 
interact with their environment* A set of hand-worked problems follows 
to illustrate and extend ideas presented in the text. Furthermore, an 
accompanying computer program, BUOY, allows the user to investigate the 
sinking rates of plankters. The module is intended for college under- 
graduates interested in aquatic biology who have some understanding of 
calculus and differential equations. 
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INTRODUCTION 

Aquatic life forms show many adaptations to the physical and 
chemical characteristics of the water In which they live. The chemical 
composition of the water is Important because the water transports minerals 
and foodstuffs dissolved in it. making them available to the Inhabitants. 
The physical properties of the water are important because they influence 
the shape, size and distribution of aquatic organisms. In this module, we 
will deal with the physical properties of still water, a portion of the more 
general subject termed fluid statics. 

Two physical properties of fluids will be of primary concern: 
density and viscosity. One effect of fluid density is manifested in terms 
of the pressure exerted on an Immersed object. The greater the depth of 
immersion, the greater is the resultant pressure on the body. Also, the 
greater the fluid density, the greater the pressure for a given depth of 
Immersion. The density of water, P^, is 775 times that of air. This creates 
J enormous hydrostatic pressures in the depths of the seas and lakes which 
creates physiological "problems" for the organisms living there. Another 
important effect of fluid density is the buoyancy force. The volume of 
fluid displaced by an- immersed object gives rise to this force, which is 
proportional to the density of the fluid. As a result, the buoyancy force 
exerted by water on an organism is 775 times that in air, the same ratio as 
the densities of water and air. Compare the structure of a tree with that 
of a large kelp plant. A tree must "use" much of its mass holding the 
photosynthetlc surfaces high to compete for light. On the other hand, 
a kelp plant being about the same density as sea water is not "concerned" 
with holding Itself up but Instead must be able to withstand the forces of 
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water currents. It is actually thicker at the top than at the bottom. 
Similarly, the greater buoyancy of water will allow most people to float, 
a feat quite impossible in the air. 

The second physical property of interest is fluid viscosity. 
The frictional resistance to movement in a fluid is due to its viscosity, 
which is caused by cohesive forces between molecules- Since the viscosity 
(y) of water is 100 times that of air, there is greater resistance to the 
povement of an object through the water as compared to movement through 
air. This increased resistance keeps plankton from sinking too rapidly 
but also retards movement of motile organisms. These animals adapt by 
assuming certain shapes for which the resistance is reduced. Thus stream- 
lined forms become more Important and more common in aquatic organisms 
than in terrestrial ones. 

In this module, then, we discuss briefly the effects of pressure 
caused by the density of water and the problems of flotation and buoyancy 
in static water. Limnologists and marine biologists seek to understand 
how pressure and flotation have influenced the characteristics and distribu- 
tion of aquatic organisms. They ask questions such as: What effects can 
high pressures have on organisms? How do organisms adapt to these high 
pressures? How do sudden changes in pressure affect organisms? In terms 
of flotation, they are interested in such things as the velocity at which 
organisms sink, adaptations that might influence the sinking rate, and 
evolutionary advantages in sinking. Many of these questions can be 
approached by examining the underlying physical principles and relating 
these to biological observations. 
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PRESSURE 

Terrestrial organisms experience a maximum pressure caused by the 
weight of the air above them of 1.01325 x 10^ dynes cm"^ (defined as 
1 atmosphere). This pressure decreases with altitude, as does the density 
because air is compressible. Aquatic organisms, however, must endure not 
only the atmospheric pressure at sea level but also an increased pressure 
caused by the weight of the column of water above them. The hydrostatic 
pressure exerted on an aquatic organism at depth h can be computed as 

P = + P^gh , (1) 

where P^ is the atmospheric pressure (at sea level given above), is the 
density of the water (1 g cm~^), g is the gravitational constant, 
980.665 cm s ^, and h is the depth below the surface (here to keep the 
units consistent, it is given in cm). The first term, P^, represents the 
atmospheric pressure at sea level and the second term represents the 
pressure caused by the mass of water above the organism at depth h. 
This hydrostatic pressure increases by approximately 1 atm for every 
10 meters of depth, exceeding 1000 atmospheres in some deep marine tranches. 
In general, the density of water is a function of both temperature and 
pressure, but these effects are small, and the valus of 1 g/cm^ is 
commonly used. 

In natural environments, fish generally experience pressures 
ranging from 0.5-500 atmospheres. Very few freshwater species encounter 
hydrostatic pressures greater than 40 atms because only 18 of the world's 
lakes have depths below 400 meters. Generally speaking, most fresh water 
fish experience even less pressure since few live below 100-200 meters. 
Marine fish, in contrast, are abundant down to 5000 meters with fewer 
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collected below this depth. Special forms have been found as far down as 
7579 meters. 

The question of how the aquatic organisms manage to adapt to these 
pressures was explored around 1891 by the French scientist Regnard. Sur- 
prisingly, he found that pressures up to 100 atmospheres had little or no 
effect on organisms and that the effects of pressures up to 1000 atm were 
often reversible — the organisms returning to normal after a short recovery 
period at normal pressures. The observed effects were functions of both 
the pressure and the time of exposure. After this pioneering work by 
Regnard, studies were neglected for 30 years until they were taken up by 
Fontaine in 1920 in France. In the 1930 's and again in the 1960's, a re-^ 
surgence in research on the effects of pressure changes occurred in Germany 
and the United States. This research determined that one of the effects 
of the increased pressure occurs at the cellular or molecular level where 
the pressure seems to cause protein denaturization. Thus, organisms living 
permanently at great depths probably have adapted by changing proteins to 
forms that are not readily denatured by the intense pressure. 

Another effect of increased pressure of biological impo?^'tance is 
related to the behavior of gases. The equation of state for a perfect 
gas is 

PV = RT 

where P is the pressure, V is the volume, T is the temperature and R is 
the universal gas constant (equal to 8.314 J g""^ mole"^ K"M- At low 
pressures, real gases tend to obey the perfect gas law. At pressures of 
less than 200 atmospheres, the error is typically less than 10% for atmos-- 

.* \ • 

pheric gases. Thus, the perfect gas law can be used to represent the 
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behavior of gases with which we will be concerned. If the temperature 
is regarded as constant, the equation (2) becomes simply 

"T PV = constant 

which is Boyles' law. The volume of gaseous spaces in organisms such as 
swim bladders in fish can change dramatically as the organisms rise or 
sink in the water. A rapid decrease in pressure can cause expansion of 
internal gases and explosion* Fish rapidly brought to the surface by a 
trawling rig can have greatly distended swim bladders. Fish that are 
forced from deep reservoirs through power turbines to shallow ponds run 
similar risks • This explosion is also experienced by divers when 
ascending too rapidly and causes damage to the lungs. Nitrogen bubble 
formation in the bloodstream, known as "the bends,'' can cause paralysis 
or even death. The increased pressure also causes changes in the quantity 
of gases dissolved in the tissues. As pressures increase, the volume of 
gases decreases, causing implosion under drastic conditions. Gases can 
come out of solution and can cause toxic effects. The most familiar of 
these effects is nitrogen narcosis. Nitrogen gas has an intoxicating 
effect on the nervous system. 

At the extreme depths of some deep lakes, the pressure may reach 
high levels (especially in comparison to terrestrial values), but on the 
whole pressure is considered a minor problem in lakes and rivers. This 
may be due not to true insensi tivity but simply to lack of research in 
that area. For example, certain aquatic insects which use spiracular 

gills for respiration have been shown to be sensitive to pressure \changes 

\ 

caused by the increase in water depth due to flooding* At high pressures 
the gills collapse or become completely wetted. Thus these sensitive species 
*3o not occur in areas where floods cause water depth changes (Hynes 1972). 
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BUOYANCY OF ORGANISMS 

As we have seen in the previous section, some fishes use swim 
bladders to maintain neutral buoyancy* Denton (1960) has described two 
other intriguing mechanisms used by marine organisms to regulate their 
depth* The cuttlefish changes its density with an osmotic pump connected 
to the cuttlebone* The animal is limited to depth above about 260 meters 
because hydrostatic pressure will overcome the osmotic potential. The 
deep-sea squid regulate density by filling the coelomic cavity with a 
fluid containing ammonium ions, a by-product o£ metabolism, which is less 
dense than sea water. This mechanism has the advantage that it is not 
influenced by changes in static pressure but it takes up a large volume 
of the animal • In contrast to these methods some fish such as sharks and 
tuna just swim all the time. This short paragraph reviews some of the 
mechanisms used by animals to stay afloat. In the rest of^^the module we 
will concentrate on smaller life forms of the marine and aquatic environ- 
ments * 

The plankton represents a community that is completely dependent 
on water as a place to live. In order for the plankton to exist in this 
unique way, they must be able to remain suspended or at least be able to 
maintain their position in the photosynthetic region long enough to grow 
and reproduce. In still water, the plankter is unable (without an external 
force) to remain truly suspended without movement but instead is actually 
sinking slowly because its density is greater than that of the water. 

Continuous residence within the euphotic zone may be neither a 
necessary nor a practical requirement for successful planktonic ^3f:^stence. 
Successful growth and survival only require that the residence time in 
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the euphotic zone permit enough photosynthesis to counterbalance respiration 
and permit some growth. The rate of growth will be influenced in part by 
the actual time spent in the euphotic zone each day. While full 'apprecia- 
tion of the total planktonic distribution must ultimately include considera- 
tion of the gross effects of water movement (waves, convection, turbulence), 
treatment of buoyancy and flotation mechanisms in static water can give 
valuable insights. 

Before proceeding we need to discuss the Reynolds number concept. 
This will help the student interpret the results from the following sections 
on terminal velocity of resistive forces on, and biological responses of 
sinking plankters. 

The Reynolds Number 
The Reynolds number. Re, is used to distinguish between two very 
different types of flow regimes — laminar and turbulent. Laminar or stream- 
line flow is characterized. by the orderly motion of fluid particles in 
essentially parallel paths. A thin dye tracer will not diffuse signifi-- 
cantly in laminar flow but will form a thin line parallel to the flow axis. 
Turbulent flow on the other hand is characterized by random motion, in which 
a dye tracer will rapidly diffuse throughout the fluid. Thus, if the 
Reynolds number is low (Re < 1.0) the flow is laminar and if it is high 
(Re > 10.0) the flow is tui:bulent (in between is a transition zone). The 
Reynolds number is defined as 
He . ^ ^ P„ 



where v is the fluid velocity, d is the characteristic dimension aach as 
the diameter of a sphere or cylinder, is the density of the water and 
y is the viscosity of the water. It was derived by non-dimensionalizing 
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che governing equations- Thus the distinction between laminar and turbulent 
flow can be made independent of size and velocity since they are included 
in the calculation of Re. Check to see that Reynolds number is a dimension- 
less quantity. 

Motion and Dynamic Equilibrium 
Consider a passive plankter (such as algae) sinking in a colurnn of 
water. The forces acting on this plankter are 1) the force of gravity acting 
downward, 2) a buoyant force equal to the weight of the fluid displaced by 
the plankter acting upward, 3) the force associatd with the acceleration 
of the plankter downward in the y direction, and 4) a force representing the 
resistance of the madium to the movement of the plankter acting upward. This 
is shown diagramatica] ly below: 

Resistive force, R- 
ma 

The buoyant and gravitational forces can be considered constant for a given 
plankter. (However, many organisms, e.g. blue-green algae, form gas vacuoles 
and can thus control buoyant forces.) Usually the gravitational force is 
greater than the uuoyant force, therefore the plankter Is sinking. The 
resistive force is not const^-int but is proportional to the velocity of the 
plankter and increases as the plankter accelerates, View'ng the spherical 
plankter as a point mass, we can s-um up the forces of the "free body diagram" 
and set them equal to zero for dynamic equilibrium (ZF=0.). Taking the down- 




Buoyant force, B 



Gravitational force, G 
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ward direction as positive we have ' 

-ma + G- B- R= 0. (A) 
To proceed with the analysis each component of equation 4 will be- discussed . 

Since the mass is equal to the density times the volume the gravita- 
tional force, G, is equal to mg = p_ V g. 

B 

The buoyant force is proportional to the fluid density, as described 

earlier.. It^^turns out that this force is just equal to the weight of the 

fluid displaced by the body which can be written p V g, where p is the 

w w 

water density, V the plankter volume, and g the gravitational constant. 

This "law of buoyancy" is known as Archimedes' principle. If the plankter 

is spherical, its volume is 4/3 Tir^, where r is the radius. Consequently, 

the gravitational force, G, and the buoyant force, B, are respectively 

A/3 "JTr^p-g and 4/3 Ttr^p g for a spherical plankter. The resistive force 
B W 

may take a variety of forms depending on the velocity at which the sphere 
is sinking. For the low vo.locities at which most plankters sink, the re- 
sistive force has been found to be directly proportional to the velocity 
and can be represented by Stokes' law, which is developed and discussed later 
in this module. Stokes' law states that the resistive force is 
R = kTTrv 



I have chosen ma as negative because the particle is decelerating. 
It is possible to assign arbitrary direction to this force but the solution 
will i-:.termine which way the object is moving. Here this choice facilitates 
a solution. 

2 

Alternatively this same equation may be derived by applying Newton's 
Second Law, which states that the acceleration of a body is directly propor- 
tional to the resultant summation of the forces acting on the body. 
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where k is a constant of proportionality equal to 67t for spherical objects, 
p is the viscosity of the fluid, and v is the velocity at which the body 
is sinking- 

Finally, recalling that the acceleration is the time rate of change 
of velocity, we can substitute each term back into equation to get equation 
3- 

(5) 



dt 3 



-™ ^ + 4 PtjS - T "^r^ p^g - 6iT |i r V = 0 



*4 

Using the volume-mass relationship again to substitute for y '^^ 



and rearranging, we have 

dv 6Tryr 



Wv 



V + g(l-T^) = 0. 

dt m 

This i.^ a differential equation of 1st order which can be solved by separa- 
tion of variables to deteimine the velocity at any time when we include the 
initial condition that at t = 0, v = 0. 
The solution is 

(6) 



V = [CgCl - -^)](1 - e"^^^) 



B 



where 



C = 



m 



6tt r u 

This solution has the general form illustrated in the following 



diagram. 



velocity 




time 



is 
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The time constant, C, which is based on the geometric and material traits of 
the plankter, characterizes the rate at which the velocity of the plankter 
approaches the limiting value- The larger the value of the time constant, 
the slower the approach and conversely as illustrated below: 



velocity 




velocity 




C large 



time 
C small 



Terminal Velocity of a Sinking Plankter 
The limiting value that is approached as the time interval gets 
larger (i.e., as t ->■<») is 



(7) 



This is called the terminal velocity, which represents the velocity at which 

the acceleration is equal to zero, i.e*, = 0* To understand this phenome- 

dt 

non, consider a body initially at rest which is suddenly released. Under the 

dv 

influence of gravity, the body begins to accelerate downward, i*e. — f 0. 

dt 

Looking at (4), one can see that this must be due to an imbalance of forces 

in this case, the resistive force is smaller than the combined gravitational 
and buoyant forces* However, as the velocity of the body increases, so does 
the resistive force, until there is zero acceleration and no change in 
momentum. The result is a constant rate of fall. 

For very small objects, the terminal velocity is approached so rapidly 
that the period of accelerational can be ignored. To demonstrate this solve 
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equation 6 for the time (t^) that it takes for the plankter to reach a 
fraction x of the terminal velocity v^. Thus 

xv^ v^(l - e ^ ) 

t /C 
e ^ = 1 - X 

t = C iln(l - x) . (8) 
c 

Consider a small plankter whose radius is r = 20 microns and density Pf 

p =1.02g cm"^. The viscosity of the water is 1 x 10"^ dyne s cm""^ . The 
B 

terminal velocity is v^ - 1.73 x 10"*^ cm s"*^. The time it takes to reach 
a velocity of 1.56 x 10^^ cm s"-^ or 90% of the terminal velocity is 2 x 10""* s. 
Thus the period of acceleration is so small that, for most computations, it 
can be ignored. 

Resistive Forces on Sinking Plankters 
As stated previously, the equation representing the resistive force 
may have a variety of forms depending on the velocity, unlike the gravita- 
tional and buoyant forces. It is known, for the examples considered here, 
that the resistive force at low velocity is a function of 1) the surface 
area of the body in contact with the water, 2) the speed of the water moving 
by the body or the speed of the body moving through the water, and 3) con- 
stants which are related to the nature of the fluid and its temperature (the 
viscosity and density of the fluid are dependent on temperature). Using 
functional notation, we write 

R = f(A, v, y, p ) 

w 

where A is the surface area of the plankter body, v is the velocity, y is 
the viscosity of the water, and p is the density of the water. 
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Assuming that our analysis of the relevant variables determining R 

has been complete, we are left with the problem of determining the exact 

functional form of f(A, v y, p ). A common tool of fluid mechanics useful 

w 

in this regard is that area of mathematics known as dimensional analysis. 
This is based upon the fact that an equation expressing a physical relation- 
ship between quantities must have the same dimensions on both sides. The 

process involves combining the independent variables (A, v, u and p ) in 

w 

Such a way that the result yields dimensions identical to those of R, Since 
the dimensions of the independent variables are in general not the sams, 
their addition or subtraction is not possible, hence the most general form 
of equation possible is a combination of products of powers of the relevant 
variables , 

In this case, A is assumed to be some power function of the radius > 
r, therefore 

R = K'r^;''p^%^ (9) 

where K' is a constant of proportionality, (Notice that independent variables 
of this function are the same ones that form the Reynolds number.) Each of 
these quantities is now replaced by its dimensional form in terms of mass, 
M; length, L; and time, T. Dimensionally , (9) becomes 

MLT-2 = (l")(l''t-^)(m"l-^")(mV^T-^) . 

Now equating the powers of each dimension on the two sies of the equation, 
we get 



ERIC 



18 



M: 1 = u + s 

L:l = n + q-3u-s 

T: -2 = -q - s 

Since there are four unknowns and three equations, we can solve simultan- 
eously for n, u, s in terms of q, and the power of the velocity. Equation 
(9) becomes 

R = K'rVp ^"V"'^ (10) 
w 

Thus, if the above assumption holds and R varies as some unknown powers of 

r, V, p and y, then a general form for the resistive force has been found, 
w 

The coefficients n, u, s were expressed in terms of q, the power of the 
velocity, because experimentation has shown that the resistive force is 
proportional to some power of the velocity. For small velocity valres , it 
is proportional to the first power; for large velocity values to the second 
power; and for intermediate velocity values some intermediate power. With 
these experimental results in mind, values for q of Ljr;..1.5, and 2 are chosen 
for discussion. 

It should be emphasized here that dimensional analysis is an experi- 
mental technique, useful in that it allows intelligent formulation of ex- 
periments leading to a maximum amount of information from a minimum of 
experimentation. The reader is referred to Vennard and Street (1975) for - 
a complete discussion. Analytical determination of R is also possible through 
a consideration of the basic equations of motion for a fluid, which is 
beyond the scope of this module. The interested reader is referred to 
Batchelor (1970). tVere an analytical solution not available, verification 
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of the usefulness of the functional form derived by dimensional analysis 
would be possible only by comparison with experiment. 

If we assume q = 1 in (10), i.e., that the resistive force is 
proportional to the first power of the velocity, then 

R = K'rvp . (11) 

This expression is the same as that derived by Stokes in 1851 by classical 
hydrodynamic (i.e. analytical) methods given earlier in the text and is 
commonly called Stokes' Law. Thus, the two methods, dimensional analysis 
and classical hydrodynamic methods, have led to the same results, indicating 
the strength of dimensional analysis. 

This relationship is important in the study of the sinking of small 
plankton because in general, for biological objects less than 0.5 mm in 
diameter, passive sinking rates fall in the range within which Stokes' law 
is valid. Bacteria, noncolonial algae, many species of colonial algae, 
nearly all protozoa, and most rotifers fall in this range, but the larger 
Crustacea and some nonspherical forms do not. Thus, if the Reynolds number 

Re = ^ ^ ^ Py < 1.0 . 

y 

then Stokes' law is valid because the flow is laminar. 

Next, if we assume that q = 3/2 in (10) then R becomes 

R - K'r W l''V^/2 . (12) 
c w 

This equation, called Allen's law, was derived in 1900 for spheres in the 
transition zone between laminar and turbulent flow (1.0 < Re < 10 .'O) . Allen's 
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law can be used to describe the resistance for large Crustacea and plankters. 

Note that the term is not the radius but computed linear dimension. 

Allen (1900) and Hutchinson (1967, p. 272) show how to compute r . 

c 

Finally, if we assume that q = 2 in (10) then 

R = K*r^v/p^ . (13) 

L W 

This relationship is credited to Newton and is valid when the flow is com- 
pletely turbulent (i.e. Re > 10). It is useful in the present context when 
considering the sedimentation of heavy materials. It also describes the 
resistive forces experienced by fish and birds as they swim or fly. In 
this connection, it is termed drag. See problem 10 for further development 
of these ideas. 

Biological Responses 

Now that the general physical principles regarding the passive 

sinking of a spherical plankter have been discussed, let^s go back and 

examine some biological responses to the problem of flotation. Examining 

Stokes' law and Allen's law, we see that the dependence of the terminal 

velocity on the excess density, (p - p ) , is very important. Indeed, if 

B w 

the plankter were able to decrease this difference to zero, i.e., make its 
body density equal to the density of the medium, the terminal velocity would 
be zero and thus the organism would not sink at all. Any reduction of the 
excess density would cause a decrease in the terminal velocity. Many 
organisms do include gases or fatty and oily materials to reduce the over- 
ail density and retard sinking, or even provide for positive buoyancy. 

Density changes on a daily basis may partially account for* diurnal 
migration* Daphnia pulex has been shown to exhibit a diurnal variation in 

51 



sinking speed with the greatest value occurring after sunrise when the 
Daphnia have finished feeding and are sinking into the deeper waters during 
the day. Many organisms exhibit density changes during their life history 
with dispersal stages or feeding stages being lighter than other stages, 
Calanus finmarchius and Calanus plimohius show gradual density decreases 
from eggs to Copepodite stages and an increase in the adult forms. It is 
interesting to note that the Copepodite feed on phytoplankton near the sur- 
face whereas the adults do not feed. Not only does density change diurnally 
and during the life history but for algal forms the density increases with 
age and disease causing older forms to sink faster out of the photosynthetic 
zone- Since the density of the medium changes with temperature, the organ-- 
isms must adapt their internal density with the season in order to maintain 
the constant terminal velocity. For example, assuming extreme temperatures 
of 0 and 25 °C, a salt water salinity of 35% and an excess density of 0.05, 
the sinking rates for a 20 micron cell calculated from Stokes' law would be 
0.023 and 0.0A5 mm s^^ at the two extremes* Similar results can be found 
for fresh water populations. It has often been assumed that cyclomorphic 
changes in some organisms or differences between the forms of tropical and 
temperate plankton may be due to adaptations to these density changes, but 
no conclusive proof exists for either confirmation or refutation. 

Empirical Results 

The above physical theory was developed to describe the sinking of 
inanimate objects such as metal spheres. Biological limnologists and 
oceanographers have shown that use of these physical principles to describe 
the sinking of biological organisms gives reasonable results. ^ 

Because of the considerable technical difficulty in obtaining re- 
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liable data, the information on the density and sinking rate of freshwater 
organisms is extremely meager. Experimental data does, however, show that 
nearly all freshwater diatoms and most phytoplanktonic organisms .except cer- 
tain species of blue-green algae are slightly denser than the medi-um and 
sink when undisturbed. Results of experiments by Fritz (1935), as reported 
in Hutchinson (1967), using Cyclotella bodanica (diameter 32-62y mean 49V) , 
Fragilaria Crotonensis (mean 68y) and Astevionella formosa (mean 68p) are 
given in the following table. 

TABLE 1. Observed and Computed Sinking Speeds 
of Diatoms in Quiet Water ^ 

Sinking Speed 

Temperature Viscosity (10~ cm s~ ) ^ 

(°C) (Poise) Observed Calculated 

Cuclotella bodanica 6 0.0147 17.0 25.0 

20 0.0101 21.0 

Astevionella formosa 6 0.01A7 8.6 12.6 

20 0.0101 11.1 

Fragilavia crotonensis 6 0.0147 59 8.6 



20 0.0101 



10.5 
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*Table from Hutchinson (1967), p. 276. 
^Calculated using Stokes' law. 

Possible reasons for the small discrepancies between observed and calculated 
sinking speeds include inadequate geometrical description of the diatoms a.d 
distrubances by small turbulent currents in the experimental appartus. 
Other work by Grim (1951) and Lund (1959), as reported in Hi- ■ chinson (1967), 
give similar ranges of values although their actual results are different. 

Most importantly. Grim (1939), as reported in Hutchinson (;i967) , con- 
cluded that living specimens of small species of Cyclotella do not sink 
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appreciably when alive but sink 7-8 m day ^ when dead or less healthy. He 
also observed that Fragilaria crotonensis sank at rates of up to 7 to 9 m 
day~^ and Synedraaaus delicatissima sank 3-20 m day"^. Thus the range for 
one species of organisms may be considerable depending on the health and 
density of the organisms, supporting previous speculation. 

The gaps in our knowledge of the zooplankton density and sinking 
rates are not as great as that of phytoplankton. For example, densities 
of many planktonic forffis of zooplankton have been computed and typically 
the excess density, p - p , lies between 0.015 and 0.025. The value for 

is W 

Gammarus pulex, given in the accompanying table, is much higher because it 
is commonly a benthic form. 



TABLE 2. Densities of Freshwater Animals 

Rotatoria (Luntz 1928) 

Euchlanis triquetral reared at 8°C 
Euchlanis tr-iquetra^ reared at 25 "C 
Euchlanis triquetra, wild, taken in summer 
Braohionus quadridentatus , tJild, taken in summer 
B. quadridentatus f. rhenanus, wild, taken in spring 
B. quadridentatus f . rhenanus, reared at 8°C 
B. quadridentatus f. rhenanus, reared at 25°C 

Crustacea (Lowndes 1938, 1942; Hamilton 1958) 
Chiroaephalus diaphanus 
Daphnia pulex 
Holopediim gibberum 
Holopediwn gibberum 
Eudiaptomus gracilis 
Gammarus pulex 
Candona Candida 

Chordata 

Gasterosteus aculeatus 
Gasterosteus aculeatus 
Bona temporaria tadpole 



1.025 



1, 
1, 
1, 
1, 
1, 
1, 



020 
025 
025 
025 
025 
020 



at 18°C 
at 18°C 
at 18°C 
at 18°C 
at 18°C 
at 18°C 
at 18°C 



1.011 at 11.3°C 
1.017 at 7.0°C 
1.0015 at 18°C 
1.014 at 18°C 
1.023 at 8.0°C 
1.066-1.088 at 20°C 
1.025 at 7.4°C 



1.003 at 10°C 
1.002 at 20°C 
1.014 at 9.6°C 



Table from Hutchinson (1967), p. 249, 
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The rate of sinking of some small zooplankton has been computed. 
For small rotifers such as Euchlanis triquetra and Brachionus quadrvdentatus ^ 
Stokes' law is valid (see section on Resistive Forces) and the results 
are given on the following page* 

TABLE 3. Observed and Computed Sinking Speeds of Rotifers 

Density (p_) Average Terminal 

^ Velocity (V^) 

Euchlanis triquetra 1.025 0,0862 cm s*"^ 

1.020 obs. 0.0658 cm s"^ 



1.020 calc.= 0.0862 x 20/25 0.0690 cm s 



-1 



- 1 



Brachionus quadridentatus 1.025 0,0495 cm s 

/. rheyianus 1.020 obs. 0.0386 cm s~^ 

1.020 caic- 0.0495 x 20/25 0.0396 cm s~^ 



Table from Hutchinson (1967), p. 278. 

Experimentation on various Daphnia species illustrates that the 
larger specimens do' indeed fall at a velocity that would be predicted by 
Allen's relationship with Daphnia pulex falling at a rate of between 0.4 
cm s*"^ and Daphnia mayna falling at a rate of approximately 0.8 cm s""^. 
Data on marine phytoplankton (especially diatom) species is abundant* 
Smayda (1970) gives a very good table of values for size and sinking rate 
observations for many diatoms and flagellates including various life stages. 
Once again, these observations agree very well with the theoretical calcu- 
lations. 

Thus the experimental work does suggest that the relationships 
developed above do hold approximately for a variety of phytoplankton and 
zooplankton species and sizes. v 
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FURTHER CONSIDERATIONS 

The development thus far has assumed a spherical plankter shape. 
In fact, many plankters are not spherical or have a gelatinous sheath to 
change the terminal velocity from that calculated above. This section 
discusses the sinking velocity for plankters exhibiting the above compli- 
cations. The effects of other factors (such as water movement) are not 
considered, but see problem 11, which considers swimming. 

Non-spherical Form 
The theory developed above was based on the case of a plankter that 
was spherical, but few aquatic organisms are spherical in form and the 
true spherical shape is conspicuously absent in phytoplankton* . One reason 
for this may be found by examining the relationships for the buoyant and 
resistive forces. The buoyant force is proportional to the volume and the 
resistive force is proportional to the surface area. Increasing surface 
area with respect to a given volume would mean increased resistance and 
slower terminal velocity. Thus, departures from the spherical form are 
important in aiding a plankter to keep from sinking. These changes from 
spherical form may involve flattening, various forms of asymmetry, forma- 
tion of colonies, and the addition of various projections such as spines, 
hairs, cilia, etc. The simplest way to deal with the effect of these changes 
in shape on the resistive force and terminal velocity is to suppose that 
the new shape can be formed by deforming a sphere of radius r^. 

The original resistance for the sphere with radius r at terminal 

s 

velocity, assuming Stokes* law (11) applies, is 

R = 67riJr v^ 
s s ts 



ERJC ^ 



-22- 



and tha resistance at terminal velocity for the new chape is 



R = 67ru(i) r V 
a "^^r s ta 



where <i> is the coefficient of form resistance which relates the change 
r 

in resistance to changes in the amount of surface area relative to the 

sphere. Since the volume and density of the body are unchanged by the 

deformation, the gravitational and buoyant forces are respectively 

D TT-^^R and p iT-|r^g. Substituting into Equation 5 with the resistance 
B 3 s^ w 3 s° 

s 

for the new body, and the velocity assumed constant, we get: 

p 4^r^g - p -^r^g - 6Try(J) r v^ =0. 
3 s^ ^w 3 s^ r s ta 

Now, solving for the terminal velocity of the new form, v^^, yields 

v " V = ^ fir (p - P )<J) ^ y" 

ta 9 ^ s ^^B ^v'^r 

Forming the ratio of the velocity of the sphere and the new body, we get 

$ = ^ . (23) 
^ ^a 

Further, we find that if a sphere can be constructed with radius d^ which 

falls at the rate v^ , then 

ta 



r 2 

where r is called the nominal radius. This latter expression in terms of 
B \ 

the radius will prove useful because the planktologist usually wants to 
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know the new terminal velocity and can measure the various linear dimensions 
to estimate d . Then, substituting in (23) we obtain 

V = 2 V 

ta r ts 
s 

Thus, a fairly convenient method exists for determining the terminal velocity 
of general nonspherical forms. The experimental results of Luntz (1928, 
1929) as reported in Hutchinson (1967) demonstrate the importance of just 
such changes from the spherical form. He showed that BvancTriorvLs quadri- 
dentatus rhenonus without spine and a density of 1.025 sank at a rate of 
0.05 cm s"^ whereas the typical form Branchionis qixadridentatus with spine 
and density also 1.025 sank at a rate of O.OAl cm s"^. 

For simple solid forms such as ellipsoids and disks, a fairly rigor- 
ous theoretical and experimental method is possible. The reader is referred 
to Hutchinson (1967, p. 260-272) for details. In general, both ellipsoids 
and disks sink slower than spherical objects of the same volume (i.e. 
will be greater than one). Somewhat surprisingly, the elongated prolate 
form is usually more efficient than the more flattened oblate form in re- 
ducing sinking velocity. It should be noted that for ellipsoids within a 
certain small class of dimensions, t{)^ < 1, so that it is not possible to 
say that a sphere will always sink faster than an equivalent ellipsoid. A 
disk shape, however, which can be regarded as the extreme form of an oblate 
spheroid, no matter what its orientation, will descend slower than an equi- 
valent sphere for Re ^ 7. For a disk one hundred times as wide as it is 
thick, the time taken to fall a given distance broadside down is 4.3 times 
that taken of a like sphere. In contrast, cylindrical organisms with 
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diameters less than 0.01 cm, oriented horizontally, will fall roughly 10 
times as fast as a sphere of equal diameter. Thus the cylindrical form 
is not efficient in reducing the terminal velocity, especially if the 
cell walls are very dense such as in diatoms. 

The latter statement brings up the important point of considering 
density differences between various geometries of identical volume which 
we might be comparing. For example, one can by reducing the diameter and 
increasing the length make a cylinder that will fall slower than a sphere 
of the same volume. However, in all likelihood, the added structural 
material would cause such a cylinder to be denser than the equivalent 
sphere, increasing the gravitational force and probably offsetting any 
gain. Similar considerations are pertinent in the case of ellipsoids, 
also. The extra surface area of an ellipsoid, as compared with an equal 
volume sphere, will probably give the former a greater density. For ex- 
ample, in the case of diatoms, part of the advantage of increased form 
resistance is obliterated by increased density since diatoms have external 
skeletons . 

In conclusion, it is important to consider the observation that 
a body moving in a liquid will tend to cake up an orientation with its 
greatest area of projection normal to the direction of motion to minimize 
the forces acting on them. A small body may have no preferred orientation. 
However, slight asymmetries in the body shape and the presence of surface 
area projections will aid the plankter in attaining the preferred orienta- 
tion with its greatrst sectional area normal to the direction of motion. 
Thus an organism will normally orientate so that the resistive foi;ce will 
be greatest for that body shape. 
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Gelatinous Sheath 

One other structural adaptation remains that warrants discussion. 
This is the gelatinous sheaths in phytoplankton and many nonmotlle planktonic 
algae such as the blue-green algae. Although the gelatinous capsules may 
increase the difficulty of zooplankton attempting to eat them, they are 
undoubtedly important in the problem of flotation. Not only does the sheath 
mean an increase in diameter and the corresponding decrease in sinking 
speed, but it also aids in reducing the organism's overall density due to 
the low density of the gelatinous sheath. 

To determine the effect of the gelatinous sheath on the resistance 
and sinking speed of the plankter, we will proceed in a manner similar to 
that of the previous section. Consider a spherical organism of radius r 
and density whose terminal velocity by Stokes' law is 

Now add a gelatinous sheath of density p to the outside of the 

s 

organism. The new radius is increased to ar where a > 1. The new volume 
is then Tir^a^ the new mass y 7Tr^[pg + (a^ - ^^^s^ density 
therefore [p + (p - p )/a^]. The new terminal velocity then becomes 

. 2 gr^a^fp - P^+-3(Pb " P^)]^-^ 
v^2 9 S W ^3 i5 s 

where [p - p + ^ (p^ - P )] represents the excess density, 
s w 3 B s 
a 

Comparing the terminal velocities thus achieved and rela-ting to 
the coefficient of form resistance (j) , we obtain I 
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a^(p - p ) + (p„ - P ) 



Assuming > 1, i.e., that the gelatinous sheath does reduce the 
sinking rate, and collecting like density terms, factoring out (a-1) and 
rearranging, we end up with 



since a > 1. Thus, in order for the sheath to be effective in reducing 
the sinking rate, the difference in density between the sheath and the 
plankton body must be at least twice as great as the density difference 
between the sheath and the medium. Indeed, if the jelly and the medium have 
the same density, the sinking speed may be reduced to any degree by merely 
increasing the thickness of the sheath. In nature, this is unlikely to 
occur since the greater thickness would interfere with the metabloic ex- 
changes of the cell with the medium. 



can reduce its sinking velocity by decreasing its density, changing shape 
or acquiring a gelatinous sheath. All or one may act to help the organism 
achieve the desired goal to remain in the photosynthetic zone until it has 
at least accomplished its objective of growth and reproduction. 

Multiple Plankters 
When two or more spheres are sinking through a liquid in the vicinity 
of each other, they cause a change in each other's sinking speed'. In small 
groups they also tend to take up specific orientations with respect to each 
other which can be predicted. Although the theory is not well developed 




> a(a + 1) > 2 



We have considered the many ways in which a single celled organism 
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for more than two spheres (even that of two" spheres is Intractable for the 
purposes of this work), some interesting observations have been made. They 
will proceed to arrange themselves in a regular polygon on the same horizontal 
plane if the cluster contains three to six spheres. Clusters of seven or 
more do not tend to form any regular grouping. The effect on the terminal 
velocity is related to the number of spheres in the vicinity, their orienta- 
tion and the distance between the spheres as compared to their diameter 
(Jayaweero et al.. 1964). Generally speaking, spheres close to each other 
decrease, the terminal velocity by an amount dependent on the angle between 

the spheres and their distance of separation, as shown below. The velocltv 
v^ 

ratio — (where v^. represents the predicted terminal velocity of the sphere 
alone and v represents the velocity with other spheres) caused by the 
mutual interaction between the sinking spheres is a linear function of 
the ratio — when the two spheres are sinking directly beside each other 
and a quadratic function of the same ratio at other angular orientations 
(Eveson. et al . . 1959). Other researchers have noted that the spheres 
tend to rotate slowly as they are falling and may separate slightly when 
thi? distance between the centers is very small (Mathew and Smith, 1960). 
Although this theory was developed for spheres under conditions that may 
no; occur for plankter, other results do suggest that the presence of other 
plankters will affect the ultimate sinking rate. Thus, when dealing with 
organisms in a natural setting, it may be expected that their terminal 
velocity will be less than that predicted by the theory given above, es- 
pecially if the density of organisms in an area is high. 
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IJHY DO ORGANISMS SINK? 

As the previous discussion and examples have illustrated, aquatic 
organisms do not have th^ same density as the medium and do indeed sink. 
Obviously continuous sinking even at very small velocities would eventually 
cause all the phytoplankton community to sink out of the photosynthetic 
zone. Field observations have shown that this is not the case. Turbulent 
movement in the photosynthetic zone no doubt plays an important role in 
resuspension. Why not just adjust the internal density to be equal to that 
of the media? In other words, why sink at all? This subject was approached 
in a very unusual way by Mu^k and Riley (1952) in what is now a classic 
paper. They not only present an interesting insight into the possible 
reasons for sinking but also provided alternative explanations for the 
presence of spines and i^^^gular shapes. 

Firstly, they suggest that an increased sinking rate would allow 
for increased absorption of nutrients from the water* They argue that if 
an organism were to remain stationary that it would soon use up the nutrients 
in its vicinity. Movement would allow renewal of nutrients in the vicinity 
of the organism and thus aid growth* By applying the physical laws of 
heat absorption of inanimate objects through diffusion and forced convec- 
tion to a theoretical analysis of phytoplankton nutrient absorption, they 
have been able to predict rates of sinking, excess density and other charac- 
teristics of the population's for waters with different nutrient levels. 
They have also argued that morphological adaptations are not necessarily 
aids to suspension but are m^jchanisms to permit twisting and vertl^cal move- 

f 

ments. They then suggestf^d that specialization can be understood as responses 
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to problems of flotation, nutrient absorption and predation. More recently 
Garvis (1976) has reformulated Munk and Riley's (1952) work to incorporate 
the physiology of the phytoplankton. It is now known that nutrient uptake 
rates are controlled by biological characteristics of the plants unlike 
the passive diffusion phenomenon that Munk and Riley modeled. Garvis 
(1976) shows that the diffusion limitations can be important for uptake 
rates at low nutrient concentrations but less so for corresponding growth 
rates. Nonetheless competitive interactions should favor small size in 
nutrient depleted waters. The inter-relationships between ecological problems 
and available morphological adaptations has led to a general theory based 
on all sizes, namely: 



Cell Size 



Large 



Intermediate 



Small 



Chief Problem 

Flotation 

Nutritional 
absorption 

Predation 



Predation 



Adaptation 

Large vacuoles — lead to low cell 
nutrient concentration and low 
cell densities 



Formation of chains which endangers 
nutrient absorption and flotation 

Chain formation: nutrient absorp- 
tion and flotation not endangered 



Parsons and Takahasi (1973), Smayda (1974), Pasciak and Garvis (1975) 
and Titman and Kilham (1976) provide more recent experimental results on 
this general topic. The interested reader will also wish to read the 
summaries provided by Hutchinson (1967 , Chapter 20) and Smayda (1970). 
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PROBLEMS 

1. a. Pressure is given in units of atmospheres and dynes crtT^ Name at 
least two other systems of units and give conversion factors. 

b. That the hydrostatic pressure increases by 1 atm for every 10 
meters at depth was given as a rule of thumb. At 10 meters what 
is the percent change in pressure as compared with the surface? 
At what depth would there be a 1 atm error using this rule of 
thumb? 

c. How will the density of water change as a function of dissolved 
solids, temperature and pressure? 

2. In the discussion section of this module, mention was made of the 

effects on fish and other aquatic organirms of sudden changes in 
the hydrostatic pressure due to changes in depth* 

l^That would be the absolute pressures and the pressure change 
experienced by a fish which lived at 25 meters and was suddenly 
brought to the surface (0 meters)? What if it were brought to the 
surface from 100 meters? 

3. Many of the higher fish, Teleostei, contain closed air bladders 

which serve as a means for reducing the overall density of the 
fish to that of the surrounding water. The air bladder regulation 
is very_„slow - so changes'^ in the gas content due to regulation will 
be ignored. Now assume that the fish in the example above had an 

air bladder which was cylindrical with diameter 1/2 cm and length 

\ 

6 cm. What is the volume at the living depth? What would be the 
volume at the surface for each case? What would be the volume 

37 



change in each case? What can you say about the relationships be- 
tween the depth changes, pressure changes and volume changes? 

The overall density of the fish consists of the density of the body 
mass of the fish and the density of the air bladder, in proportion 
to the volume occupied by each component. Since the volume of the 
air bladder changes with pressure and depth, the overall density of 
the fish changes also. 

Imagine a salmon at a depth of 100 meters. The volume of the 
air bladder is 5 cm^ and the remaining body mass volume is 95 cm^. 
The density of the gases in the air bladder is .00125 gm/cm"^ and 
that of the remaining body mass is 1.053 gm/cm"^. Confirm this. 

If the salmon is now raised 50 meters, what would happen to the 
overall density? VHiat would this mean to the fish? I^hat if instead 
the fish were lowered 50 meters? 

In the above examples the fish had no means of controlling the 
volume change. What if instead of the salmon we placed a diver at 
100 meters? If his lung volume is 500 cm^ at 100 meters, what would 
be the volume at the surface? What is the volume change? How can 
he deal with this volume change? 

What disadvantages are there to swim bladders? What kinds of fish 
would you expect not to have them? 

Assume that you are studying a population of small spherica^ plankter 
that occur in the ocean. You collect: a sample of them and measure 
the radius and density. The average values you obtain are .001 cm 
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radius and density p = 1.001. Assuming the standard values for p 
and y what would be the terminal velocity if Stokes law for the re- 
sistance can be assumed to hold? If the plankter continued to sink 
at this rate for a full day, how far would it sink? If the radius 
were 10 times as big, what effect would it have on v^ and how far 
would it fall in a day? Would we expect this plankter to keep sink- 
ing under natural conditions? What might influence the distance 
that the plankter sinks? 

Using equations (7) and (8) compute v^ and t^ for 90% of v^ for the 

given values of p and r. 

B 

= 1.0001 1.001 1.0001 
r = .005 .005 .05 



Sketch curves of v and t vs. r and p„ - p . 

t c B w 

With r constant what does a 10 fold increase in p„ - p do to v 

B w t 

and what does the corresponding change in p do to t ? With p 

Be B 

constant what does a 10 fold increase in r do to v and t ? 

t c 



For the water temperatures of 1, 15, 25 and 40 ^C, what is the 
density and the viscosity? Which is more important for determining 
the sinking velocity? At what temperature would a plankter sink 
fastest? 

In equations 11 through 13 the resistance force is derived using 
dimensional analysis. Assuming that the acceleration is zero for 
equation 4 in the text, derive the terminal velocity for Stokes', 
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Allen's and Newton's resistance terms « Assume a spherical plankter 

Also leave the solution in terms of r for Allen's form. 

c 

Not all plankters are passive organisms. Many possess medns of 
swimming which may help in maintaining or at least regaining posi- 
tion. Such organisms as the small crustacea may be quite powerful 
swimmers whereas others have only limited success with cilia or 
flagella. For example, species of dino-f lagellates have been shown 
to have swimming speeds ranging from 0.05-3 cm min~^. 

VJrite down a force balance equation to include the swimming term 
S, using Stokes law. If the acceleration is zero, what is the 
terminal velocity. 
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PROBLEM SOLUTIONS 

1. a. Table of Conversions 
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From Halliday, D. and R. Resnick, 1962. Physics, part II- Appendix 
G, p. 44. 

1. b. A : 10 meters the change in pressure is 0.980665 x 10^ dynes/cm""^ of 

96.78% of the 1 atmosphere estimated. The depth where there is a 

1 atm error using this rule is calculated as 

.0322 atm error _ 1 atm error 
10 xa X m 

X = 310.56 m 

1. c. Density Increases with increasing pressure and Increasing dissolved 
solids. Distilled water has a minimum density at about A ^C.^ 
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Pressure 

2) P = + p^gh 



1.01325 X 10" dynes/cm^ + 1 g cm"^ x 980.665 cm x 25 



X 10^ cm = 



P23 = 3.A649 10^ dynes/cm^ = 3.42 atm 



= 1.01325 X 10^ + 0 = 1.01325 x 10^ dynes/cm^ = 1 atm 
P25 - Pq " 2.452 x 10^ dynes/cm^ = 2.42 atm 

P^QQ= 1.01325 X 10^ + 9.8 X 10^ = 10.82 x 10^ dynes/cm' = 10.68 atm 
^100 ~ ■^O ^ dynes/cm^ = 9.68 atm 



3) TTr^ = 7T(.25)^ = .196 cm^ , irr^Jl = .196(6) = 1.18 cm^ 

0 



25 

PV=C P'V =P*V V= — =— V 

' 25 25 0 ^0 • 0 P« 25 



4.04 cm^ = -1.18 , AV = 2.86 cm' 

1.18 = 12.60 cm^ 
AV = 11.42 cm^ 



Pq "^100 ''0 1 



ap ap, 

^—^ = ^rrr^ = .846 .'. AP=.846AV 

^^25 ^^100 



^ = ^=.10 AP=.10AD 
^^25 ^°100 
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^5 ^air ^air +Vxp -Vx p 

bladder bladder ^"'^^ ^"^^^ "t^- 

-5 X .00125 + 95 X 1.053 = 100 x p 

^total 

.00625+100.04 . , nn / 3 
Ptotal 100 = ^'^^ Sn/cm 

P^Q = 1.01325 X 10^ + 4.903 x 10^ >= 5.9166 x IQ^ dynes/co^ 

«= 5.839 atm 

^100 ^100 " ^50 ^50 

^50 = '5.139'^ = 5-^^ ^otal = = 

9.14(.00125)+100. 04 „^ , 3 
Ptotal 104714 = Sm/cm . 
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The fish is now lighter than the water and will rise. 
'l50 



P^^„ = 1.01325 X 10^ + 14.71 X 10^ = 15.723 x 10° dynes/cm' 



= 15.52 atm 



P V = P V 
^100 100 150 150 



V,^_ = — fl^ll^ = 3.44 cffi^, ^ = 95 + 3.44 = 98.44 cm^ 
150 15.52 total 

3.44(.0012 5)+100.04 , ,3 
^total 98744 = ^'^^ S^/^^ • 

The fish is heavier than the water and will siiit^ 

5) V P = V P 

0 0 100 100 

V3 - lO^MCim . 53A0.0 cm' 

AV = A8A0 cm^ 

The diver must decrease the volume in his lungs by exhaling deeply as he rises 

3 

so that his lung volume wil] remain at 500 cm . This is in contrast to the 

fish which cannot easily increase or decrease the air bladder ^vblume voluntarily. 
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6) A fish must change the volume of its swim bladder to achieve neutral 
buoyancy. This would be a disadvantage if the animal were changing 
depths quickly as does a tuna. Fishes inhabiting shallow waters might 
not need them. 

Buoyancy 

g = 980.665 cm/s^ y = .01 g/cm-s 

d = .001 cm p_ = 1.001 g/cm^ 
(Pn-P„) = -001 g/cm^ 

J5 W 

v^ = 21.793 d^ 

= 2.1793 X 10"^ cra/s d = .001 

Distance travelled in 1 day = v^ x // of seconds in 1 day 

= 2.1793 X 10-^(8. e/ixlO**) = 1.883 cm 
v^ - 2.1793 X 10"^ cm/s d = 0.01 

Distance = 2.17^3 x 10"'' (8 . 64x10"*) 
= 188,3 cm. 

Thus a 10 times incriiase in d causes a 100 times distance increase in 
both the terminal velocity, v^, and t:he distance. 

Under natural conditions, turbulent water conditions would cause the 
plankton to be carried to the upper surface again. 

The plankton could also change its density to change the sinking rate, v^ 
and the distance it falls. 

8) v^ - 2,792.556 d^(p„-p ) t = 5l.l7p d^ 

t B W C D 

a) p„ = 1.0001 p -p = .0001 d = 0.005 

= 5.45 lO-^cm/s 
t = 1.28 10-^s 
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1.001 P---P = 0.001 d •= 0.005 

B w 

5.45 X lO-** cm/s 
1.28 X 10"^ s 

1.0001 P^-P„ = -0001 d = .05 

5.45 X 10"^ cm/s 
1.28 X 10-^ s 

100 fold increase in for 10 fold increase 
in d 

"a 




10 fold increase in for 10 fold increase in 



100 fold increase in for 10 fold increase in d 




10 fold Increase in for 10 fold increase in p^. 



45 



-41- 



9) 



XJater temperature (**C) 


1 


15 


25 


40 


Density (g cm"') 


0.9996 


0.9993 


0.996"^ 


0.9923 


Viscosity (g cm~^s~^) 


0.01724 


0.01138 


0.00890 


0.00650 



10) 



Viscosity is more important than density in determining sinking rates. 
A plankter will sink faster in warmer water. 

Using Stokes' law and substituting into equation A with acceleration zero 
we have 

mg - P Vg - K'rv^p = 0 . 
w t 

Assuming that the plankter is spherical, this becomes 



- 67Trv^lJ = 0 ♦ 



Solving for v^ , we get 



^t = 



g(l - — ) = 9 gr (Pb-P„)U 



which is identical to the solution found in the module, 

Likewise for Allen law, we get 

mg - p Vg - kr V ^ V 3^ p U"^"^ - 0 
^w ^ c t t ^w ^ 



V, = 



1/2 r g'^^ (P^-PJ^^ P„-^' Vi^ 

C 13 w 



W 



where 



= r - 2/5 d 



and 



9y^ 



- = 2gp^(p^.p^) 



Similarly for Newton's law, we find that 



-1 



w 



Notice that in all cases the terminal velocity is dependent on the 



term P„~P called the excess density. 
B w 
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11) ™ TT = mg — ii— = 67rryv - s 

dt 



t 9y p„ 6TTry 
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Synbol 



a 
C 
D 
6 



h 
k 



n 
P 



q 

R 
R 

Re 



APPENDIX I. SYMBOLS, UNITS AND DIMENSIONS 
Quantity Units Dimensions 



area 

constant 

constant 

distance between sinking spheres 

force of gravity 

gravitational constant 

height 
coefficient 

coefficient 

mass 

exponent 
pressure 

atmospheric pressure 
exponent 
resistance 
gas constant 

Reynolds number 



cm 

dimensionless 
seconds 

cm 

Dynes 
980. cm s 



_2 



cm 

dinensionless 
dimensionless 
g 

dimensionless 
dynes cm"^ 
dynes cm"^ 

dimensionless 
dynes 

8.314 J g-' 
mole"' r' 

dimensionless 



MLT- 
L 



2 



ML-'r^ 



S.t. Equivalent 



m s 



kg 



N m' 



_2 



N in' 



-2 



U> 

I 
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Symbol 



Quantity 



Units 



Dimensions 



S.I. Equivalent 



r 



resistance of new shape 
resistance of sphere 
radius 

nominal radius 

critical radius (Allen's law) 

mks of sphere 
swiniing force 
temperature 
time 

time to reach 
fraction x exponent 
volume 
velocity 

terminal velocity 

terminal velocity at new form 

terminal velocity of sphere 
terminal velocity of 1 
terminal velocity of 2 



S 
T 
t 
t 

c 

u 
V 



ta 



ts 



'tl 
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dynes 

dynes 

cm 

cm 

cm 

cm 

dynes 

seconds 
seconds 
dimensionless 

3 

cm 

cm s' 
cm s" 

cm s" 

cm s' 
cm s" 
cm s" 



1.-2 



L 
L 
L 

L 

MLT' 

e 

T 
T 



1-2 



ir 

LT" 
LT" 

LT" 

ir 

LT" 



m 8 
t s" 
n s" 

m s' 
ra s" 
m s" 



I 

^ 

I 
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fraction 



density 



density of the body 
density of gelatinous sheath 
density of water 
viscosity 

coefficient of form resistance 



dlmensionless 

.3 



g cr 



g cm 



■3 



g'cm"^ 



g cm 



g cm"' s"' 



dlmensionless 



ML- 



ML 



ML 



MLT^ 



kg m' 
kg m' 



-3 



-3 



kg m' 



kg m"' s* 
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COMPUTER PROBLEMS 



Program BUOY is described in detail in the section titled "User's 
Guide to Program BUOY." It contains a line printer graphics subroutine 
named PRNT3D (Gales, 1978) which displays all of the output generated 
by the program. In each exercise the user must input certain parameter 
values according to the rules and restrictions listed in the user's 
guide. 

1* Plankter can alter its sinking rate by adopting different strategies 
of size and density. Consider, for example, a hypothetical plankter 
of radius = 0.01 cm. How does the theory, explained in the accompany- 
ing text, predict that the terminal velocity varies with density (use 
density values between 1.09 and 1.2)? In particular, is there any 
difference in the functional relationship for the regions where Stokes' 
law or Allen's law are valid? Why might these differences occur? 

2. Further exploring the strategies of planktonic organisms in regulating 
their sinking rates, consider a hypothetical plankter of density 1.01 
g/cm^. How does theory predict that the terminal velocity varies with 
radius (use radius values between .01 cm and .05 cm)? Are there any 
differences in the functional relationships for the regions when Stokes' 
law and Allen's law are valid? 

3. Combining the results from questions (1) and (2), find a set of radius 
and density values that results in a terminal velocity of 8.7 x 10"^ 
cm/s. Verify that these values are correct. Find other combinations 
of radius and density values which yield the same terminal velocity. 
These results suggest that an infinite number of combinations of radius 
and density will give a given terminal velocity. In nature, only a few 
of these combinations are actually seen. !-Jhy? 

4. Many hydrodynamics texts plot the graphical representation of radius or 
density versus terminal velocity on log paper. This is done typically 
because the range of values is very large an d^'SQsiffr detail is lost. Usin 
the log plot option and the values used for questions (1) and (2), com- 
pare the current plots with those obtained previously. 

5. In the previous example, the water temperature and, consequently, the 
viscosity of the water were considered constant. In nature, the temper- 
ature and viscosity change with the season of the year and with -depth. 
Consider the hypothetical plankter; find a set of radius and density 
values such that in the summer the plankter doesn^t sink any faster than 
2.0 X 10"^ cm/s. What effect does this combination of values have on 
the sinking rate during the winter? Now choose a set of values that 
results in a winter sinking rate of approximately 2.0 x 10"^ cm/s. l^Jhat 
effect does the combination of values have on the summer sinking rate? 
What, strategy is generally adapted by aquatic plankters? 




USER^S GUIDE FOR PROGRAM BUOY 



Identi f Ication 

BUOY -A Program Which Displays the Effects of Plankter Radius, 
Plankter Density, and Water Temperature on the Terminal 
VelociLy of a Sinking Plankter. 



Authors -Christina £• Covan and Mark Bailey, Center for Quantitative 
Science in Forestry, Fisheries ^ and Wildlife, University 
of Washington, Seattle. May 1979. 

Purpose 

Program BUOY is the computer supplement to the module 
"Pressure and Buoyancy in Aquatic Ecosystems" by Christina E. Cowan. 
The program is designed to Let the student investigate the effects 
of plankter radius, plankter density, and water temperature (acting 
only through water viscosity and not water density) on the sinking 
velocity of the organism. Two different water resistance formula- 
tions, associated with the names of Stokes and Allen, are used depend- 
ing on the value of the Reynold's number. 

Operation 

Several explicit and implicit facriors control the operation of 
program BUOY by way of the formulation of the resistance function and 
the selection and format of the output. 

The resistance function depends on the value of the Reynt)ld*s 
number* which in turn is a function of user-input. Stokers law is used 
if the Reynold's number is less than 0.5; Allen's law is used if it 

er|c * 55 



-48- 

exceeds 0.5 • 

A list of user-assigned options named OPLIST controls the 
selection of output. Each option in the list may take on the values 
1 through 10 and the list may contain up to 10 unique elements. The 
'Options are divided into four classes, as follows: 

Option 1: Perform calculations for a single (radius, 

density, temperature) triple and print out the associated 
Reynold's number, resistance, and terminal velocity. 
Options 2, 3, 4: Hold density and temperature constant 
and plot radius versus a) Reynold's niimber (option 2), 
b) terminal velocity (option 3), and/or c) resistance 
(option A). 

Options 5, 6, 7: Hold radius and temperature constant and 

plot density versus Reynold's number (5), terminal velocity 
(6), an:d/ "^r re^yiBtancB (7)^ 
Options 83 9, 10: Hold radius and density constant and plot 
temperatture versus Reynold*s number (8), terminal velocity 
(9), and/or resistance (100 . 
For example, OPLIST = 1,3,5,7 invokes options 1,3,5, and 7 in a single 
run. Variables held constaiit for a set of plots are specified by a 
single value, whereas variables which form a plot axis are specified 
by two values and indicate a value range. For example, in options 
5, 6, or 7, the input 

RADIUS = 0.07, 
TEMP = 10, 

HENSTY = 1.09, /TQ/ 1.2 
indicate that RADIUS (0.07) and TEMP (10) are to be held constant. 
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whereas DENST^ is to form a continuous x-axis which varies from 1.09 
to 1.2, 

BUOY calls a printer plot package named PRNT3D to generate 
.the plots for options 2 through 10. The plot formats (e.g. size, shape, 
scale (log or linear) etc.) are controlled by a large number of user- 
assignable parameters which are described in the User's Guide for 
PRNT3D (Gales, 1978). These parameters may be included along with 
the above mentioned input for BUOY; however, BUOY itself assigns de- 
fault values which are adequate for most purposes (one common exception: 
setting NY = -45 instead of + 45 generates plots with a logarithmic 
y-axis scale) . 
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Prpgram Organization 

The program is organized around the following flow chart; 



Output 
appropriate 
error 
messages 



7^ 



Yes 



T 



< NODFLT = TRUE?~> 



Yes 



No 



Read default values for 
all variables from the 
built in default file 



Read in the next 
user-^supplted data set 



< FINIS = TRUE? > ■ 



Yes 



> 


No 


Check for errors in 
the input set iust read 







Terminate 
program 



^ Errors found? y 



No 



4^ 



(^Is OPLIST exhaustedY") 



Yes 



No 



< OPLIST(I) = 1? > 



Yes 



Calculate values 
for a single point 
and print out 



No 



Calculate all points 
in a plot, and flag 
any er. or conditions 



Yes 



EKLC 



< Errors? > 



No 



Write out x,y coordinates 
plus plot titles. Call the 
printer plot routine QQI R3D 
which generates the plot 



< I = 1 + 1 > 
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Input 

All input is handled by a format free input package (Gales 
and Anderson, 1978) which permits a user to assign values to variables 
by a "name=value" convention. Not all variables need be explicitly 
assigned by the user, however, as unassigned variables automatically 
assume default values. The input consists of any number of data sets, 
each of which is terminated by a dollar sign ($). Comments in an input 
set are enclosed within slashes (/). Each data set may generate a 

number of separate printer plots. 

The input for BUOY is divided into three classes: a) variables 

having biological significance or related to program options: NPOINT, 

OPLIST, RADIUS, DENSTY, and TEMP; b) variables which control certain 

program operations, such as program termination or the handling of 

default input: IPRINT, ECHO, NODFLT, and FINIS; and c) variables which 

control the printer plots (default values are in parentheses): XMIN (0), 

XMAX (0), YMIN (0), YMAX (0), ZMIN (0), ZMAX (9), XRICH (0), YRICH (0), 

DFAULT (0), OVPRNT (.F.), AVE (.F.), INT2D C-^O, NX (60), NY (A5) , and 

ZMAP (0,1,2,3,4,5,6,7,8,9). The variables in the first two classes are 

explained in the following INPUT TABLE, whereas the printer plot variables 

are explained in rhe user's guide for PRNT3D (Gales, 197S). 
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INPUT TABLE 



^Name 


Type and 
Dimensions 


Range 
Limits 


Description 


NPOINT 


Integer 


[2,99] 


The 
for 
is : 


number of points calculated 
each curve. Defaul-t value 
NPOINT « 12 


'OPLIST 

f 

i 

f 

i 


IntegerdO] 


[1,10] 


A list of up to 10 unique options 
which control the output. The 
options are as follows: 
1: Perform calculations for a 

single (RADIUS, DENSTY, TEMP) 
triple and print out the 
associated Reynold's number, 
resistance, and terminal 
velocity 








2: 


Hold (DENSTY, TEMP) constant 
ana p±ou xviUJXUo vs • x\eyno±a s 
number 


i 

! 

! 
1 






3: 


Hold (DENSTY, TEMP) constant 

SLTLQ pXOu KAUJLUo VS • LcXulHiaX 

velocity 


1 






4: 


Hold (DENSTY, TEMP) constant 
and plot RADIUS vs. resistance 


i 




i 

1 


5: 


Hold (RADIUS, TEMP) constant 
and plot DENSTY vs. Reynold's 
number 


i 
1 
r 






6: 


Hold (RADIUS, TEMP) constant 
and plot DENSTY vs. terminal 
velocity 


1 






7: 


Hold (RADIUS, TEMP) constant 
and plot DENSTY vs. resistance 








8: 


Hold (RADIUS, DENSTY) constant 
and plot TEl'iP vs. Reynold *s 
number 


• 






9: 


Hold (RADIUS, DENSTY) constant 
and plot TEMP vs. terminal 
velocity 








10: 


Hold (RADIUS, DENSTY) constant 
and plot TEMP vs. resistance 








If OPLIST(I) = 0 for some I, then no 
action is taken for the T-th option* 
Default values are: OPLIST = 
1,0,0,0,0,0,0,0,0,0 
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Input Table Cont. 



Name 


[ Type and 
Dimensions 


Range 
Limits 


Description 


RADIUS 

! 
1 

1 


Real (2) 


[0.001,1.0] 


The radius of a plankter, in centi- 
meters. For options 1,5,6,7,8,9 
and 10, RADIUS is specified by a 
single number (or the first number 
only in a pair); for options 2,3 
and 4, it is specified by a pair 
of numbers which indicate a range 
of values (e.g. RADIUS =0.1, /TO/ 
0.3). Default values are: RADIUS 
= 0.01, /TO/ 1.0. 


1 DENSTY 
1 

i 

1 
) 
! 
i 
I 

1 
1 

i 
i 


Real (2) 

1 

I 

1 


[1.001,2.0] 


The density of a plankter, relative 
to water (= 1.000). For options 
1,2, 3, A, 8, 9 and 10, DENSTY is speci- 
fied by a single number (or the 
first number only in a pair) ; for 
options 5 , 6 and 7 , it is specified 
by a pair of numbers which indicate 
a range of values. Default values 
are: DENSTY = 1.09, /TO/ 2.0 


i 

TEMP 

1 

! 1 

I 

i i 


Real (2) 

1 


[0,30] 


The water temperature, in degrees 
centigrade. For options 1,2,3,4, 
5,6 and 7, it is specified by a 
single number (or the first number 
only In a pair); for options 8,9 
and 10 it is specified by a pair of 
numbers which indicate a range of 
values. Default values are: TEMF 
= 0.0, /TO/ 30 


i IPRINT i 

; i 
f 

' 1 
I 

! i 

! i 


Logical 

i 


. F . , # T . 


A Logical value which causes the 
current values for all input var- 
iables (default as well as current 
user input) to be printed. Default 
value is: IrRINT = .F. 


I i 
! ECHO ! 

i i 

i 


Logical 


• F . , . T . 


A Logical value which causes the 
user's input to be echoed if ECHO 
= .T,, or suppresses echoing if 
ECHO = .F. Default value is: 
ECHO - .T. 


NODFLT 


Logical 


.F., .T. 


A Logical value* which suppresses 
the input of def -^ult va'lues if 
NODFLT = .T. Default value is: 
NODFLT = .F. 
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Input Table Cont. 





Type and 


Range 




Name 


Dimensions 


Limits 


Description 


FINIS 


Logical 


• F • > • T , 


A Logical value which causes program 








termination if and only if FINIS 








= .T. Default value is^: FINIS = .F. 
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The last four variables deserve special mention. 

1. The logical variable IPRINT controls the output of all input variables 
which are currently in effect (default values as well as those speci- 
fied in the current input set). Setting IPRINT=TRUE (or T or ,T.) 
displays the current values for all Input variables; setting IPRINT= 
FALSE (or F or .F.) suppresses the display. 

2. The logical variable ECHO controls the echoing of the input cards. 
Setting ECHO=TRUE causes the subsequent input set to the echoed; 
setting ECHO=FALSE suppresses the echo for the subsequent input set. 

3. The logical variable NODFLT can be used to inhibit the automatic 
assignment of default values to input variables. If NODFLT is set 
TRUE in the current input set, then the current input set is assigned 
default values as usual » but all s ubsequent input sets merely accumu- 
late more input values. In effect, the input values which exist after 
the i-th input set is read, become the default values for the (i+l)-th 
input set. The standard default values may then be restored by setting 
NODFLT=FALSE, but, again, the effects of this change are delayed until 
the next input set is read. To a limited extent, NODFLT permits a 
user to set up his own default values and can be very useful for exe- 
cuting a number of input sets which differ only in a few parameters* 
Consider the following example in which a user wishes to produce a 
number of radius versus terminal velocity plots with different values 
for temperature: 

/INPUT SET 1: THE FOLLOWING VALUES BECOME THE/ : 
/DEFACTO DEFAULTS FOR ALL SUBSEQUENT INPUT SETS:/ • \' 

NODFLT = TRUE, ??POINT = 15, OPLIST = 3, 

RADIUS = 0.001, /TO/ 0.5, TEMP = 0, DENSTY = 1.3, $ 
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/INPUT SET 2: UP THE TEMPERATURE/ 

TEMP = 2, $ 
/INPUT SET 3: UP IT AGAIN/ 

TEMP = 4, $ 
/INPUT SET 4: AND AGAIN/ 

TEMP = 6, $ 
/INPUT SET 5: FRY IT TO A CRISP/ 

TEMP = 30, $ 
/INPUT SET 6: STOP IT NOW/ 

FINIS = TRUE, $ 

The logical variable FINIS controls program termination. The user 
must add the card: 

FINIS = TRUE, $ 

as the very last input set. If FINIS is not set, the program will 
terminate abnormally . 

Output 

BUOY produces a small table of values for option 1, and a maximum 
of nine plots for the full set of options 2 through 10, in a single run* 
Each plot contains a set of titles at the top of the plot, x and y axis 
annotation, x and y axis numeric labels, and a plot legend which speci- 
fies scale factors for the numeric labels. The parts of a plot are 
best explained by an example. Consider the first plot in the set of 
sample runs* The plot title consists of the four lines: 
HYDRAULIC PROPERTIES OF PLANKTON 
RANGE OF X (.0100, .0400), RANGE OF Y (.0138, .8805) 
DENSITY = 1.0100 
TEMPERATURE = 0 
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The X- and y-axis annotations are: 



REYNOLDS NUMBER 



(y-axis) 



RADIUS (CM) 



(x-axis) 



The X and y numeric labels, along with the first line of the plot legend, 



point in the extreme upper right corner (designated with a "1" character) 
is at coordinates (x = 4.000, y = 8.805). Since the scale factor for x 
is E - 02 (i.e. 10 ) and the scale factor for y is E - 01 (i.e. lO"" ) , 
the true coordinates are (x = 4 x 10 , y = 8.805 x 10 ) « (x = 0,04, 
y = 0.8805). Note that the ^y-axis is plotted on a logarithmic scale. 
The last two lines in the plot legend show the number of points mapped 
to each of 10 z-axis' levels (this information is more useful for programs 
which generate a number of different curves on a single page, each one 
identified by a single z-axis level). In this case, a large number (> 99) 
were mapped to z level 0 and show up as blanks, 20 were mapped to z level 
1 and show up as I's, and none were mapped to any other z level • 

Restrictions 

The user must restrict all input variables according to the range 
limits listed In the input table, in order to avoid unrealistic physical 
values. 

Error Messages 

There are four types of errors which may occur when attempting to 
execute program BUOY: 

1. Syntax errors in the user's input 

2. Range check errors \^ 

3. Calculation errors 

4. Plot parameter errors 



enable the user to interpret the plot numerically. For example, the 
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For type 1 and 2 errors, the program flags the error, skips the calcula- 
tions and plotting, and then reads the next input set. For type 3 errors, 
the program continues with the calculation and plotting after o? tputting 
an error message. For type 4 errors, the program suppresses plotting, 
outputs an error message, and reads the next data set. For a complete 
description of type 1 and type 6 error messages and actions, refer to 
the user's guides for the format free input package and printer plot 
package, respectively. 

Range check errors occur if the range limits set for input variables 
are exceeded. These messages are of the form: 

ERROR NO. X IN BUOY 

yy • • • y = dd • • • d 
yy **'y = dd***d 
ONE OR MORE OF THE ABOVE 
VARIABLES IS OUT OF RANGE 
where x is an error number from 1 to 4, yy^^y is a variable name (NPOXNT, 
RADIUS, DENSTY, TEMP, OPLIST) , and dd-.-d is a value. 

The calculation error results when illegal combinations of input 
values escape the range checks. It is of the form: 

ERROR NO. 5 IN BUOY 

THE REYNOLD'S NUMBER dd— -d 
IS GREATER THAN 10. 

Sample Runs 

The annotated listing starting on the next page illustrates the 
control cards and input cards for several sample runs. Each input set 
is terminated by a $ and generates 0 or more printer plots. ^ 

66 

ERIC 



59- 



XBUO'Y* T15# P2. 
ACCQUNT#9XM0eC00» 

comment! ♦♦♦*♦♦♦♦♦»♦♦*♦**♦♦♦*♦♦***♦♦♦*♦♦»**************♦****♦****♦**** 

comment! ♦ THE FIRST CARD SPECIFIES THE JOB NAME, -XBUOY-, THE ♦ 
COMMENT. ♦ CENTRAL PROCESSOR TIME, -T20- (IN SECONDS), AND THE JOB ♦ 
COMMENT. ♦ PRIORITY, -P2-. THE SECOND CARD SPECIFIES THE BUDGET AND ♦ 

COMMEN.T. 

ATTACH^BBUOYjID-BBUOY . 

attach, bff, id-bff. 
attach,bpr3d, id-bpr3d . 

comment! ♦♦♦♦♦♦♦*♦♦*+****♦***♦***♦♦**♦♦*♦♦♦♦♦*♦♦***♦****♦**♦****♦♦♦♦♦♦ 

comment! ♦ THE ABOVE CARDS ATTACH THE PROGRAM BUOY, -BBUOY-, AND THE * 
COMMENT. ♦ SUPPORT ROUTINES FREE FORM INPUT, -BFF-, AND PRNTBDj * 
COMMENT. ♦ -BPR3D-. THEY ARE ALL IN BINARY FORM. 

COMMENT. ********* *******^************* ********************* ********** 

COMMENT . 
LOAD, BBUOY. 
L0A0*3FF. 
LaAD,BPR3D. 
EXECUTE, BUOY. 
C OMMENT . 

COMMENT. ♦♦♦♦♦♦♦♦♦*♦*♦♦♦***♦***♦♦*♦*♦********************************* 

COMMENT. * THE ABOVE CARDS LOAD THE PROGRAM AND SUPPORT ROUTINES, * 
COMMENT. * AND P/iSS CONTROL TO BUOY FOR EXECUTION. 

COMMENT. ♦♦♦♦♦*♦♦*♦♦**♦*****♦**♦♦**♦♦♦♦♦*♦*♦****♦♦***♦***♦************ 

COMMENT . 
EOR 

THE FOLLOWING DEFAULT VALUES ARE ASSUMED UNLESS OVERRIDDEN BY 
INPUT VALUES 

OPLIST • 1 

NPQINT • 12 

DENSTY - 1,09,2.0 

RADIUS ■ 0. 001, l.-O 

TEMP ■ 0.0, 30.0 

ECHO - .T. 
FINIS - .F. 

NOOFLT « .F. 

IP«INT ■ .F. 

♦♦♦♦♦***♦♦*♦♦***♦♦** ************RUN 1* 

DPLIST-1, R ADIUS-0.007, DENSTY-l.Ol, S 



♦♦♦♦♦*♦♦*♦♦*♦**♦♦♦**♦♦♦********* RUN Z*********** 



********* 



OPLIST-2, 3/ A/ NP0INT-2D, RADIUS-O.Ol, /TO/ 0.0<r, NY- -'t5, 
OENSTY-l.Ol, S 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦**♦** ♦♦♦♦♦♦♦♦♦♦♦♦RUN 3^ ♦♦♦♦*♦♦♦♦♦♦♦*♦♦♦♦***** ******** 

ERIC S7 ' ' 



0PLlST-5#6#7# NPDINT-15* DENSTY-1.09# /TO/ 1.2* t 

/ 
/ 

/ ***r^*r^*****i^fi\JH i^* *************** *************** 

/ 

0PLIST»8*9# 10> NP0INT«12# DENSTY-1.01# RADIUS-0.0^* S 

/ 
/ 

/♦♦♦♦»♦«♦♦♦♦♦♦♦♦♦♦»♦* ♦•♦♦♦♦♦♦STOP PROGRAM* 
/ 

F-INIS«.T.?S 

*EDR - 



\ 
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JGRAy -BUOY- READY FOR INPUT 

' / THE FOLLOWING DEFAULT VALUES ARE ASSUMED UNLESS OVERRIDDEN BY 
INPUT VALUES 

OPLIST - 1 
NPOINT ■ 12 
DENSTY ■ 1,09,2.0 
RADIUS - 0.001/ 1. 0 
TEMP « 0.0/3C.0 
ECHO • .T, 
FINIS ■ .F., 
NOOFLT - .F. 
IPRINT • .F . 

^«4> *«<l>*4>**4 ♦4' + + + + * + + RUN ♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦***♦♦*♦»♦♦ 

DPLIST-l* RADIUS-0.007/ OENSTY-l.Olt S 

US .70000E-02 

;ITY elOlOOE+Ol 
'Ew.AiURE - 0 

iOlOS NUM3ER .*t71P0E-02 

STANCF .1^069E-0^ 

IINAL VELOCITY .5996^E-02 



J- 
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Pft!D6RA.J) -BUOY- READY FOR INPUT ^ 

I ' I 

/ / 
OPLlS7-2j3#<t^ NPDINT-20> PADIUS-0.01, /TO/ 0.0<i/ NY- -<»5> 
DENSTY-l.Olf $ 
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HYDRAULIC PROPERTIES OF PLANKTON 
RANGE OF X ( .0100> * ,'*{^§|oo'' ^ ^ •0138» .8805 ) 

TEMPERATURE** 0 



1.000 1.^58 1.966 2.<t75 



2.983 



3.^92 



6.805 



5.A89 



3.^22 



2.133 



1.330 



517 



,322 



,201 



^.000 

— X 
i 



13BY1 

X — 
. 0 00 



1.^53 



1.966 



2.A75 2.983 
RADIUS (CM) 



. — X > X 

3.«t92 ^.000 



SCALE FACTORS ■ X-AXIS» E-02 

0(-9), 1.000(20)* 
Z5-Z9* 5.0C0( 0)* 6.000( 0)* 



Y-AXIS(LOG) « 
2.000( 0)* 
7.000( 0)f 



E-01 
3. 000 { 
8. 000 ( 



Z-AXIS» E+00 
0), <i.000 ( 0) 
0), 9.00 0( 0) 
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RAMGE OF X 



HYDRAULIC PROPERTIES OF PLANKTON 
( .0100* .0<.00)> RANGE OF Y { .0122* 
DENSITY ■ 1.0100 
TEMPERATURE « 0 



.1665) 



1.000 
X — 

1.665' 



1.238 



l.'.SS 
X — 



1.966 
— X— • 



2.<i75 



.920 



. 65AY 



3 7H 



.281 



.209 



ire 

t K J J 



.122 Yl 

V 

1 .000 



2.983 
X 



3.A92 
X 



4.000 
• — X 

1 



— X — 
1 .4b6 



— X 

J .966 



— — yt, — — — — — -——Y, — ~' 
2.<t75 2.98 3 

RADIUS (CM) 



SCALE FACTORS - X-AXIS I E-02 
ZO-2'i- 0(-9 ), 1.000 (20), 

25-Z9- 5.000( 0)# 6.00c; 0)# 



Y-AXIS (LOG) I 
2.000( 0)t 
7.000( 0)» 



E-01 
3. 000 ( 
8.000 ( 



X X 

3.A92 4.000 



Z-AXISI E+00 
0 ) , 4 .000 ( 0 ) 
0), 9.000( 0) 



?2 
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I HYDRAULIC PROPERTIES OF PLANKTON 

. RANGE OF X ( .0100. .0A00)» RANGE OF Y ( .0000» .0A87) 

DENSITY ■ 1.0100 
I TEMPERATURE - 0 



1.000 
X— 

A.37AY 



2.180 



1.A58 
X— 



1.966 
— X — 



2.^75 



2.983 
X— 



3 .A92 
• — X — 



.975 



,^36 



195 



,0b7 



.039 



.017 



.00^- 



.0G4Y1 

X— 
1.000 



4.0C 

' — X 
1 



• — X — 
1.^58 



— X— 
1.966 



z'lld ",963' 
RADIUS (CM) 



■ — >; 

3.^.92 



■ — X 
'♦.OO 



mi 



FACTOR^, 
Z5-Z'J- 5.000( 



■ X-AXISJ 
•9), 1.0 

0), 6.000 ( 



E-02 
000 (20), 
0), 



Y-AXIS (LOG ) I 
2.000( 0), 
7.000{ 0)f 



E-C2 
3.000 ( 
8.000 ( 



F + QO 
A.000( 0) 



Z-AXIS I 
0 ), 

0)f 9.000( 



0) 
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PfiOGRAIJ -BUOY- READY FOR INPUT ^ 

/ ' / 

/^^ ^ 4 *************** **fi,{iH 2* *********** w ****************** / 

/ I 
OPL ^ST-5*6# 7» NP0INT»15> DENSTY-1.09, /TO/ 1.2, 1 
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HYDRAULIC PROPERTIES OF PLANKTON 
RANGE OF X { 1.0900> 1.2000)* RANGE OF Y ( .X238> .2752) 

RADIUS ■ .0100 
TEMPERATURE ■ 0 



1.090 
X— 

2,752Y 



2,580 



1.107 

— ^-X 



1. 125 



l.l^A 



1 .163 
X- — 



i.iei 

— X— 



R 2.236 

£ . 

Y 

N 

0 

L 2.06A 

0 



N 1,<3?? 

U 

M 

B 

E 

R 1.720 



1. 5A8 



1.376 



1.233 



1.200 

• — X 
1 



1 
X — 
1.090 



, — X 

1.107 



— X 

1.125 



— X _.-.x— — — - — — X—- — — X 

1.1^^ . 1.163 l.iei 1.200 

DENSITY (GM/CM**3) 



SCALE FACTORS ■ X-AXISi E+00 
ZO-2^- Oi-9), 1.000(15)/ 

25-Z9- 5.000( 0)* 6.000( 0), 



Y-AXISi E-01 
2,000( 0), 

7.000( 0), 



Z-AXISt E+00 
3.000{ 0)* ^.000( 0) 
8.000( 0)* 9.000( 0) 



ERIC 
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HYDRAULIC PROPERTIES OF PLANKTON 
RANGE OF X ( 1.0900/ 1.2000)/ RANGE OF Y (. .1101* .2443) 

RADIUS ■ .0100 
TEMPERATURE • 0 



1.090 
X— 



2.295 



1.107 
• — X — 



1. 125 
— X— 



1. l<i^ 
— X 



2. 1<»2 



1.989 



1. ?36 



1.683 



1.530 



1.377 



1.22A 



1. lOlYl 

X — 
1 .090 



1.163 
X— 



1 .IBl 

X— . 



1 .200 



■X 



■ — X- — 
1.107 



— X — ■ 
1. 125 



■ X——— —————— X — —————— X — — — — X 

1.1^^ 1.163 1.181 1.200 

DENSITY (GM/CM**3) 



SCALE FACTORS - X-AXIS« E+00 
ZO-Z^» 0{-9)/ 1.000(15)/ 

Z5-Z9- 5.000( 0)/ 6.000( 0)/ 



Y-AXIS I 
2.000( 
7,000( 



E-01 

0)/ 

0>/ 



Z-AXISt E+00 
3.000( 0)/ 4.000( 0) 

e.oooi 0)/ 9.ooo( 0) 
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HYDRAULIC PROPERTIES OF PLANKTON 
RANGE OF X ( 1.0900* 1.2000)* RANGE OF Y ( . 000<»» . .0008 ) 

RADIUS - .0100 
TEMPERATURE - 0 



1.090 



1.107 
V ... 



1 .125 



1. 1^^» 



1 .163 



1.181 



8 .203Y 
I 
I 



7.691Y 
I 
I 
I 
1 

7. 173Y 
I 
I 
I 
I 

6.665 Y 

T 



1.200 



6.152 



5.6A0 



5.127 



^ .61^ 



^1,102 



3.691 



1 
X— 
1 .090 



— X- — 
1.107 



— X — 
1.125 



SCALE FACTORS ■ X-A XI S t E+00 
Z0-ZA« 0(-9), 1.000(15), 

Z5-29- 5.000( 0), 6.000( 0), 



X — — — X— " — X — — )( 

l.l^A 1.163 1.181 1.200 

DENSITY (GM/CM**3) 

Y-AXISI E-OA Z-AXISI E+00 
2.000( 0)f 3.000{ 0)* A.000( 0) 
7.000( 0), 8.000{ 0), 9.00O( 0) 
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PROGRAf^ -BUDY- READY FQP INPUT ^ 

/ - / 

/ ^ If ************************ ******RUH *,*******************************/ 

I / 
0PLIST"e#9^ 10* NP0INT-12> DENSTY"1.01» RADIUS»0.0^» t 
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HYDRAULIC PROPERTIES OF PLANKTON 
RAi^GE OF X ( 0*30.0000)/ RANGE OF Y < .8805* <r.3187) 

RADIUS • .0^00 
DENSITY - I. 0100 



319 



3.928Y 



3.537Y 



3. 1^*7 



2.755 



2.365 



1.97AY 



1.193 



966 
• X 



1.983 
• — X— 



361 



3^J00 
1 



'1 



-X- — 
.^58 



.-X X— X- ~—X— — — — — X 

.966 1.^75 1.983 2.49? 3.0C0 

TEMPERATURE (CELSIUS) 



SCALE FACTORS ■ X-AXISt E+01 Y-AXISi E+00 2-AXISi E+00 

0(-9)/ 1.000(12)* 2.00G( 0)» 3.000 ( 0 )t '♦.000( 0) 
Z5-29- 5.000( 0), 6.000( 0)* 7.000( 0)* 8.000( 0)* 9.000( 0) 
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HYDRAULIC PROPERTIES OF PLANKTON 
RANGE OF X ( 0*30.0000)f RANGE OF Y ( .1665» 

RADIUS • .0^00 
DENSITY ■ 1.0100 



.966 



l.<f75 



1.98 3 

■ — — 



.3HS0I 



2.^92 



3. 160Y 



2.990Y 



2. 820Y 



N 2 .650Y 
A 



2. A80 



2.311 



M 2.1^1 

I 



1.971 



1.801 



1.665Y1 



3.000 
• — X 
1 



-X-- 



■ X X X--- 

966 1.^75 1.983 

TEMPERATURE (CELSIUS) 



SCALE FACTORS « X-AXISi E+01 Y-AXIS« E-01 
Z0-2^t» 0(-9), 1.000 (12)f 2.000 { 0)» 

Z5-29« 5.000( 0)/ 6.000( 0)* 7.000( 0»* 



Z-AXIS S 
3.000 ( 0 ).> 
8.000 ( 0 )v 



. — x--« 
2.^92 



E + 00 

^.000 ( 

9w000 C 



• — X 
3.000 



0) 
0) 
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RANGE OF X 



HYDRAULIC PROPERTIES OF PLANKTON 
( 0»30,0O00)» RANGE OF Y ( .0^87* 

RADIUS ■ ,0<tQO 
DENSITY ■ 1.0100 



1. 896Y 
I 

I 



.^58 
-X 



.966 
-X — 



1.475 

, — X 



1.983 
X-- 



.1896) 



X 



3.000 
• — X 
1 



1.736Y 
I 



- 1 



I • 
1.576Y 
I 
I 
I 
I 

1.A16Y 
I 
I 
I 
I 

1.256Y 
I 
I 
I 
I 

1 .096Y 
I 
I 
I 
I 

.936Y 
I 
I 
I 
I 

.776Y 
I 
I 
I 

I 

.615Y 
I 
I 
I 

.<»87Y1 
X- 



-X~- 
.458 



X — x — ————— — X — — ■- X 

.966 1.^75 1.983 2.492 3.000 

TEMPERATURE JCELSIUS) 



SCALE 

Z0-2A- 

Z5-29r^ 



FACTORS - X-AXTSi E+01 Y-AXiSl E-01 
0(-9)* 1.000(12)> 2.000( 0)f 
5.000( 0), 6.000( 0)^ 7.000( 0)/ 



Z-AXISJ E+00 
3.000{ 0)f 4.000{ 0) 
6.000{ 0), 9.000( 0) 
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/PROGRAf^ -BUOY- READY FDR IMPUT 

/ / 

*********** **********'^4^^$yQP p^QQ^/^n********' 

I 

FINIS • .T., J 

PROGRAM -BUOY- TERMINATED 
EOP 
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